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具有 行为 效应 的 含 AGV 柔性 车 间 调 度 研究 ` 
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Jj 要 : 为 了 研究 行为 效应 对 含 AGV 有 柔性 车 间 调 度 问题 的 影响 ， 构 建 了 考虑 学 习 效 应 以 及 恶化 效应 的 以 完工 时 间 最 
短 为 目标 的 数学 模型 。 针 对 该 模型 的 特点 , 采用 基于 工序 、 机 床 和 AGYV 分 配 的 编码 方式 , 并 提出 具有 混沌 局 部 搜索 策 
略 的 鲸鱼 算法 求解 。 数 值 算 例 结果 验证 了 算法 的 有 效 性 和 可 行 性 。 实 验 对 比 了 学 习 效 应 、 恶 化 效应 、 学 习 - 恶 化 效应 、 
不 考虑 效应 四 种 情况 下 加 工时 间 的 差异 。 同 时 研究 了 加 工 任务 批量 、AGV 数量 对 加 工时 间 的 影响 。 通 过 实验 证 明了 
AGV 资源 约束 和 各 效应 共同 影响 了 单位 批量 完工 时 间 。 
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Research on flexible Job-Shop scheduling problem with AGVs constraints and 
behavioral effects 
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(School of Business, University of Shanghai for Science & Technology, Shanghai 200093, China) 


Abstract: In order to study the effect of behavioral effects on flexible job-shop scheduling problem with AGV, this paper 
constructed a mathematical model with the goal of shortest completion time considering learning effect and deterioration effect. 
According to the characteristics of the model, this paper adopted a coding method based on processes, machine and AGV 
allocation, and proposed a whale algorithm with chaotic local search strategy. The results of numerical examples verify the 
effectiveness and feasibility of the algorithm. The experiment compared the processing time differences between the learning 
effect, the deterioration effect, the learning-deterioration effect, and the absence of the effect. At the same time, this paper studied 


the influence of the number of processing tasks and the number of AGVs on the processing time. Through experiments, it is 


proved that the AGV resource constraints and various effects jointly affect the unit batch completion time. 
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问题 更 加 符合 实际 地 生产 制造 。 当 机 器 重复 加 工 同一 零件 ， 随 

着 加 工 数 量 的 增加 ， 加 工 单个 零件 的 时 间 缩 短 ， 即 产生 学 习 效 

AGV Cautomated guided vehicle) 是 一 种 可 以 通过 智能 计算 应 ; 当 零 件 开 工时 间 变 晚 ， 使 加 工 单个 零件 时 间 变 长 ， 即 恶化 
H 


机 控制 系统 进行 控制 的 物料 搬运 移动 式 机 器 个 。 它 在 生产 制 效应。 文献 [3] 提 出 具有 恶化 和 学 习 效应 及 资源 分 配 的 排序 模型， 
造 领域 应 用 最 广泛 ， 因 此 以 AGV 搬运 约束 的 柔性 车 间 调 度 为 ” 且 将 模型 一 般 化 为 户 =P maxUf (0,5) eet - p (其 中 已 为 工件 
背景 去 进一步 深入 研究 具有 现实 意义 。 j 的 可 用 资源 ，e 成 为 恶化 率 ，max{ 刀 7 学习 效应 函数 )。 文 
在 实际 生产 过 程 中 ， 工 件 加 工时 间 会 受 行为 效应 的 影响 。 献 [4] 研 究 的 是 考虑 学 习 恶 化 效应 的 无 等 待 流水 调度 学 习 恶化 
三 义 的 行为 效应 是 通过 研究 人 的 行为 规律 、 机 器 加 工 规律 ， 找 ”效应 模型 ， 利 用 局 部 搜索 方法 求解 最 小 总 延迟 。 文 献 [5] 考 虑 了 
到 提高 生产 效率 的 途径 ， 包 括 学 习 效 应 、 恶 化 效应 等 。 本 文 研 ”同时 具有 学 习 恶化 效应 的 柔性 作业 车 间 调 度 问题 ， 并 用 改进 的 
究 工 人 具有 学 习 效 应 与 工件 具有 恶化 效应 时 对 加 工时 间 的 影响 ，“ 遗传 算法 结合 可 变 区 域 搜索 进行 求解 。 考 虑 到 参数 设置 会 影响 
则 在 帮助 企业 更 合理 地 安排 生产 加 工 计划 。 文献 [2] 考 虑 学 习 恶 ” 算法 性 能 ， 采 用 了 田口 的 鲁 棒 设计 方法 来 定义 最 佳 参 数值 。 文 
化 因素 并 且 机 器 具有 可 用 性 限 钥 i 


1， 实际 加 工时 间 为 献 [6] 考 虑 学 习 效应 以 及 设置 时 间 的 柔性 车 间 调 度 问题 ， 提 出 遗 
P,-(p,*btyr* (2 为 工件 恶化 率 ，4 为 学 习 因 子 ) 使 得 单机 调度 ， 传 算法 结合 变量 领域 搜索 和 从 代 局 部 搜索 ， 并 采用 自 适 应 策略 
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来 调整 参数 。 文 献 [7] 考 虑 具 


学 习 恶 化 效应 的 平行 机 调度 问题 ， 


提出 了 一 种 基于 蚁 群 优化 算法 求解 。 文 献 [8] 考 虑 了 同时 具有 学 
习 恶 化 效应 的 柔性 作业 车 间 调 度 问题 ， 并 用 改进 的 遗传 算法 结 


合 邻 域 搜索 进行 求解 。 


笔者 在 查阅 文献 的 过 程 发 现 : 具有 行为 效应 的 含 AGV 4X 


学 习 恶 化 效应 的 调度 问题 ， 


生 车 间 调 度 研究 极 少 。 目 前 学 者 研究 考虑 学 习 遗 忘 效应 ， 或 者 


但 研究 者 这 往往 取 一 个 国定 的 学 习 


率 、 遗 忘 率 、 恶 化 率 求解 分 析 与 不 含 效应 的 调度 问题 对 比 然后 
就 结束 了 研究 。 往 往 忽略 分 析 不 同 的 学 习 率 ， 不 同 的 遗忘 率 、 


不 同 的 恶化 率 的 影响 差异 怕 


E。 本 文 不 仅 考 虑 了 上 述 不 同 的 效应 


因子 对 加 工时 间 的 影响 ,还 分 析 了 任务 批量 与 AGV 的 数量 对 含 
AGV 柔性 车 间 调 度 的 影响 。 


1 含有 AGYV 的 FJSP 描述 与 建 模 


1.1 问题 描述 
含有 AGYV 的 柔性 车 间 


调度 问题 考虑 n 工件 在 m 台 设 备 上 


加 工 并 由 AGV 负责 工件 在 设备 间 的 运输 的 过 程 ， 且 每 台 设备 


都 由 一 个 工人 操作 。 每 个 零 


件 都 有 i 道 不 同 的 工序 m<m, 每 个 


工序 可 在 多 台 设 备 上 加 工 ,有 T 台 AGV, H r<m+2。 有 如 下 假 


设 ; 


a)AGV 的 搬运 路 线 是 固定 的 ,AGYV 搬运 不 延 时 ,不 同 AGV 


Us ACTU. 
b) AGV 工作 内 容 是 : 


从 上 一 道 工 序 加 工 设备 的 输出 缓存 


区 立即 搬运 至 下 一 道 工序 的 输入 缓存 区 。 


c) 零 时 刻 ， 所 有 设备 和 机 器 都 可 用 
d) 工 件 从 加 工 设 备 的 输入 缓存 区 到 开始 加 工 的 时 间 是 瞬时 


完成 ， 消 耗 时 间 忽略 不 计 


e) 每 台 设 备 一 次 只 能 加 工 一 个 零件 , 每 个 零件 在 设备 上 加 


工时 不 中 断 。 
Ð 设备 加 工 顺序 约束 ， 
则 。 


加 工 设 备 按 先 到 先 加 工 (FCFS) 规 


g) 不 同 工 件 先后 无 约束 , 同 种 工件 必须 先 完成 上 一 道 工 序 ， 


然后 再 进行 下 一 道 工序 。 


ho 所 有 待 加 工 的 工件 放 在 装载 台 M0， 加 工 好 的 工件 放 在 


EAG Mm-1 处 。 
1.2 ”模型 建立 
1.2.1 学 习 效 应 模型 


学 习 效 应 是 指 在 实际 生产 中 ， 工 人 随 着 加 工件 数 的 增加 ， 


操作 逐渐 熟练 使 得 加 工 工 从 


的 时 间 缩 短 。 又 因为 学 习 效果 因 人 


而 异 ， 使 得 每 个 人 加 工 相同 工件 时 间 不 同 。 因 此 很 有 必要 研究 


学 习 效应 对 加 工时 间 的 影响 。 
参数 的 取 值 可 能 使 实验 结果 出 现 较 大 偏差 ， 然 而 目前 尚未 


有 针对 学 习 效应 参数 取 值 的 研究 。 因 此 本 文采 用 简洁 的 DeJong 
学 习 效应 模型 四。 模型 如 下 : 


P, - PIM & 0- M)r^] (1) 


其 中 : P, 


效应 时 工件 的 加 工时 间 , r 表示 工件 排 在 第 r 位 加 工 。 表示 学 


a=- (2) 


Td 


表示 实际 加 工时 间 , PB 表示 理论 加 工时 间 , 即 没有 学 习 
oeh 


习 因 子 。 


1 表示 学 习 率 。M 表示 不 可 压缩 因子 : 当 M=0 时 ， 即 


Biskup 模型 u0， 当 M=1 时 ， 则 表示 工件 加 工时 不 考虑 学 习 效 
应 的 情况 。 


工件 实际 生产 时 间 


图 1 学 习 效应 示意 图 


参考 文献 [11, 12] 可 知 ， 员 工学 习 效果 受 员工 的 初始 能 力 、 
员工 的 能 力 上 限 、 员 工 的 学 习 能 力 和 重复 作业 次 数 四 个 方面 的 
因素 影响 。 起 初 员工 到 岗位 作业 时 拥有 不 尽 相 等 的 初始 能 力 ， 


而 且 员 工 的 初始 能 力 和 员工 工作 经 验 、 是 否 有 岗 前 培训 有 关 。 
某 一 岗位 员工 的 初始 能 力 用 员工 加 工 工件 的 时 间 与 标准 加 工 工 


时 的 比值 表示 , 初始 能 力 越 强 ，? 越 小 ; 工作 越 复杂 ，2 越 大 。 
构建 员工 初始 能 力矩 阵 pw (其 中 py G-12, mj =1,2,...,n) 


表示 员 


i 加 工 工件 j 的 初始 能 力 ， 一 人 独立 操作 一 台 加 工 设 


备 ， 操 作 加 工 设 备 共 灵台 ， 工 件 ” 共 件 ): 


in 
Pu Ph Pi 
initial initial initial 
LI Pai P» P; 
initial initial initial 
Pm Pm Pon 


学 习 能 力 与 员工 的 观察 力 、 记 忆 力 、 理 解 能 力 、 注 意 力 等 


员工 本 身 有 关 。 员 工学 习 能 力 用 学 习 率 1 (Learning Rate? 


表示 ， 学 习 能 力 越 强 ， 那 么 学 习 率 1 小 ;工作 越 简 单 ，! 越 大 。 


构建 员工 学 习 率 矩阵 L( 其 中 (i=42…,1m;j=1,2,…,n;) 表示 员工 


i 加 


随 着 重复 操作 次 数 变 多 ， 工 件 加 工时 间 趋 于 稳定 ， 即 员工 


能 力 达 至 


件 j 的 学 习 率 ): 


mn 


| 上限。 员工 能 力 上 限 与 工作 本 身 难 度 、 员 工本 身 有 关 ， 


工作 越 难 , 能 力 上 限 越 低 ; 员工 学 习 能 力 越 弱 , 能 力 上 限 越 低 。 
其 中 pi" G-152, m j 212,0) 表示 员工 i 加工 工件 j 的 最 大 能 
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力 。 
final final finat] 
11 12 ssi ln 
final final final 
pio" 一 | ja 22 tee 2n 
final final final 


ml m2 s m | 


重复 作业 次 数 越 多 ， 员 工学 习 效 果 越 好 ， 用 r 表示 员工 第 
r 次 加 工 该 工件 。 根据 学 习 效应 的 四 大 影响 因素 修正 DeJong 学 
习 模 型 得 : 


Ba esu --M)r*] [E pP" P, 
ür 


ü 
final final 
Pi P; Fp < py P 

lgl; 
E 
lg2 


reN* 


(3) 


为 了 验证 员工 的 初始 能 力 、 员 工 的 能 力 上 限 、 员 工 的 学 习 
能 力 和 重复 作业 次 数 是 影响 员工 技能 变化 的 重要 因素 ， 假 定 有 
两 个 员工 ， 分 别 加 工 同 一 批 工 件 同一 道 工 序 : PY =0.9， 
1,208, p^"207 p" =-08，1 站 =09，Phm -075, M-05. 
两 人 技能 曲线 变化 图 如 下 图 所 示 。 由 图 可 得 ， 员 工 2 初始 技能 
强 于 员工 1, 但 是 学 习 能 力 弱 于 员工 1， 当 重复 加 工 5 次 时 ， 员 


Tur 


工 2 的 技能 弱 于 员工 1， 并 在 t=6 时 ， 两 人 技能 保持 不 变 。 
^g — -uri| 


0.65 上 


重复 加 工 次 数 r 


图 2 员工 技能 曲线 图 


1.2.2 恶化 效应 模型 

实际 生产 加 工 中 ， 当 工件 的 开工 时 间 足 够 晚 或 者 位 置 足够 
后 面 时 ， 工 件 的 加 工时 间 不 可 能 无 限 大 。 然 而 ， 在 我 们 能 看 到 
量 的 关于 恶化 效应 研究 成 果 中 ， 都 忽视 对 这 一 方面 的 研究 。 
因此 本 文 在 一 般 线性 恶化 模型 的 基础 上 进行 模型 改进 。 一 般 线 
ESLIRE, P =P tbt 其 中 三 为 基本 加 工时 间 ， 为 该 工 
件 开 始 加 工时 间 ，0< 己 <1 为 恶化 因子 且 只 与 工件 有 关 。 


P? = P, +b [max (t 207^ ,0) — max (t 2£7^,0)] 
| Uu ij J J J (4) 


B=[b b, .. b] 


n 


O 
3. 
d 


AREA, 


等 : 具有 行为 效应 的 


Pij*(tmax-tmin)bj 


[ 件 实际 生产 时 间 


tmin 工件 开 工时 刻  tmax 


S 


3 恶化 效应 示意 


1.2.3 学 习 - 恶 化 效应 模型 


学 习 效 应 主要 是 工人 通过 技能 的 熟练 ， 使 加 工时 间 缩 短 ， 
然而 恶化 效应 是 指 工件 开工 时 间 延 后 导致 加 工时 间 变 长 。 加 工 
时 间 长 短 受 两 者 共同 作用 ， 因 此 有 必要 综合 考虑 学 习 效 应 、 恶 


化 效应 对 工件 加 工 的 影响 ， 即 学 习 - 恶 化 效应 。 

本 文采 用 Cheng (2017) 提出 学 习 - 恶 化 模型 0 : 
F,-(*ba)r? ” ， 结 合 本 文 改 进 的 学 习 效应 模型 得 改进 学 习 - 恶 
化 模型 : 


P, N p" (P, +b [max{t -e ,0} max (t = 1? ,0} DIM 十 a Nu M)r” ] 
H "B, > pP" (P, tbt) 
P, = pP" (P, +b [max(t — t?" ,0} — max (t —(*,0)]) 
HP, < pP" (P, +b,t,) (5) 
lgl; 
[^ 二 一 二 
lg2 


reN* 


2 ”鲸鱼 算法 设计 


2.1 算法 编码 和 解码 


本 文 用 一 个 座 头 鲸 位 置 代表 一 个 候选 调度 方案 ， 其 中 维 数 
X 3L, H Tr (transport sequence) 向 量 、M (machine) 向 量 、 


R (robot) 向 量 构成 ， 且 这 三 个 向 量 每 一 维 上 是 相互 对 应 的 。 则 


Lenin Cn 为 工件 数 , 为 第 i 个 工件 的 总 工序 数 ), H Tr 


向 量 分 量 小 于 等 于 工件 总 数 、M 向 量 分 量 小 于 等 于 加 工 设备 
CRRA HRE) 总 数 、R 向 量 分 量 小 于 等 于 AGV 的 数 
Eo X 1 为 一 个 编码 示例 。 

表 1 粒子 编码 向 量 矩 阵 


ap 


地 


粒子 维 数 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
Tr 向 量 1 2.3 4232 14 1 3 2 3 4 1 3 
M 向 量 1 2 3415231 5 4 6 3 6 6 6 
RH 1 2 12 12 121 2 1 2 1 2 1 2 


按 鲸鱼 优化 算法 更 新 位 置 公式 计算 ， 不 可 避免 出 现 非 整 数 
值 。 又 因为 在 FJSP 中 工序 有 先后 顺序 和 工序 、 加 工 设备 、AGV 
必须 为 整数 。 因 此 按 如 下 规则 调整 迭代 结 

a) Tr 向 量 。 由 算法 公式 可 知 , 分 量 不 可 避免 会 出 现 小 数 。 
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针对 FJSP 中 工序 有 先后 顺序 约束 和 工序 为 整数 的 特点 ， 采 用 


下 图 方法 修正 迭代 结果 。 首 先 将 迭代 后 的 Tr 向 量 按 数 值 大 小 从 


小 到 大 到 排序 后 得 到 新 的 向 量 ， 然 后 将 新 向 量 的 值 所 对 应 的 初 


始 值 重 新 排列 ， 最 后 得 到 修 


可 行 工 序 加 工 方案 。 


正 后 的 向 量 序列 ， 即 得 到 一 个 新 的 


初始 Tr 向 明 1 2 E334 2 
ERG 21 7 [285 


E 3 4 3 1 3 


2.9 03 29 0.8 -0.1 0.9 0.8 23 -05 87 79 23 


从 小 到 人 24 7 27 24 29 03 了 9 798.01 09 08,427 -05 8&7 79 23 
-0.5 -0.1 03 08 08 09 21 21 23 HSEQREI29 29 7 79 87 
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1 1 4 s R2 2 2 1 4 


图 4 


Tr 向 量 调整 


b) M 向 量 。 由 算法 公式 可 知 , 分 量 不 可 避免 会 出 现 小 数 。 
针对 FJSP 中 加 工 设备 为 整数 的 特点 (加 工 设备 满足 全 部 柔性 )， 
采用 下 图 方法 修正 迭代 结果 。 首 先 将 迭代 后 的 M 向 量 向 上 取 


整 ， 因 为 加 工 设 备 数量 的 限定 ， 然 后 按 边界 取 值 ， 即 将 所 有 超 
出 上 界 〈 加 工 设 备 总 数 ) 的 取 值 为 上 界 , 低 于 下 界 的 取 值 为 1。 


若 加 工 设 备 为 部 分 柔性 ， 首 


先 获 取 已 修正 的 Tr 向 量 , 得 到 加 工 


工序 顺序 。 然后, 将 迭代 后 


加 工 设备 集 的 加 工 设备 编号 距离 最 近 的 数 更 新 为 新 的 M 向 量 。 


M 向 量 的 分 量 与 对 应 工序 可 选 加 工 


初始 MI 量 | 2 3 4 1 5 2 3 | 5. 4 6 3 6 6 6 


迭代 后 32 8 45 24 48 0. 


向 上 了 到 整 
| 按 边 界 到 位 


图 5 


3 29 0.8 -01 56 09 27 23 87 79 23 


M 向 量 调整 


c)R 向量。 调整 规则 同 


HRH 1 2 3 4 1 5 2 3 1 5 4 6 3 6 6 6 


iMm 24 8 56 21 34 0. 


M 向 量 ( 加 工 设 备 满足 全 部 柔性 )。 


3 29 0.8 -0.1 56 09 -3.0 23 97 -98 67 


BERE 3 8 6 3E 
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BURRE 2 2 2 2 


图 6 
22 种群 初始 化 


R 向 量 调整 


根据 上 文 编码 方式 ， 第 一 步 随机 生成 搬运 序列 ， 第 二 步 采 


用 贪 禁 法 分 配 机 器 , 第 三 步 启 发 式 生 成 R 向 量 , 并 借鉴 文献 [15] 
启发 式 调整 的 思想 ， 调 整 R 向 量 ， 使 空 载 时 间 减 少 。 


采用 贪 梦 法 分 配 机 器 : 


a) 随机 生成 的 Tr 向 量 可 得 到 


ds 


加 工 工序 顺序 。 


b) 逐一 计算 工序 的 最 日 
配给 对 应 加 工 设备 加 工 。 


完工 时 间 , 则 将 最 早 完工 的 工序 分 


c) 若 多 台 加 工 设备 的 完成 工时 间 相 等 , 那么 判断 该 工序 加 


入 是 否 需 要 等 待 , 那么 分 配 等 待 


时 间 短 的 加 工 设备 加 工 该 工件 ; 


若 等 待 时 间 相 同 ， 随 机 分 配 一 台 加 工 设 备 加 工 该 工件 。 


启发 式 生 成 、 调 整 R 向量: 


a) 生成 R 向 量 : 为 每 个 搬运 序列 的 搬运 工序 启发 式 均匀 


分 配 AGV 搬运 ， 得 到 R 向 量 


FH. 


boys R 向 量 : dE AGV 搬运 工件 到 达 加 工 加 工 设备 后 ， 
检查 该 加 工 设备 缓存 区 是 否 有 工件 等 待 搬运 。 如 果 有 ， 则 该 


AGV 紧 接着 搬运 缓存 区 的 


件 ， 即 将 工件 对 应 的 搬运 工序 调 


整 到 该 搬运 操作 之 后 , 从 而 减少 该 AGV 的 空 载 时 间 ; 如 果 无 ， 
则 执行 搬运 序列 的 下 一 个 搬运 操作 。 根 据 R 向 量 ， 将 搬运 序列 


从 左 至 右 逐 一 检查 ， 从 而 调整 R 向 量 ， 最 终 使 R 向 量 的 AGV 
空 载 时 间 最 短 。 
2.8 位置 更 新 策略 
鲸鱼 优化 算法 06( whale optimization algorithm,WOA) 是 由 
亚 学 者 Mirjalili 等 人 于 2016 年 时 提出 的 一 种 新 型 群体 
智能 优化 算法 , 其 源 于 对 自然 界 中 鲸鱼 群体 捕食 行为 的 模拟 ; 通 
过 鲸鱼 群体 搜索 、 包 围 和 追捕 攻击 猎物 等 一 系列 过 程 实现 搜索 
优化 的 目的 。 虽 然 WOA 已 被 证 明 在 求解 精度 和 收敛 速度 方 范 
均 优 PSO、DE、 引 力 搜 索 算 法 (gravitational search algorithm, 
GSA) 和 GWO 算法 等 , 但 其 依然 存在 全 局 搜索 与 局 部 搜索 难以 
ii. BERN Ul, Tp WOA 算法 在 这 些 问题 
上 的 改进 研究 国内 外 仍然 尚未 展开 。 
WOA 算法 思想 可 知 ， 座 头 鲸 的 位 置 是 “ 游 走 砚 食 ”“ 包 
所 收缩 ”“ 螺 线 捕食 ”三 种 行为 分 互相 影响 作用 的 结果 。 
a) 游 走 员 食 。 座 头 鲜 依 靠 种 群 随机 个 体位 置 来 导航 搜索 猫 
物 ， 其 空间 位 置 更 新 如 下 公式 所 示 。 
=X ma -DA 


E 
zt 
Es 


su 


H 


1+1 


D-|cC- X.. ,- Xl (6) 


其 中 : X, Xa 分 别 表 示 座 头 鲸 第 t. tH 次 迭代 的 位 置 。A、C 
为 系数 向 量 且 Ce[0,2] ，Ae[-2,2] ， 当 4<1 时 ， 表 明 随 机 个 体 
为 猎物 所 处 位 置 ， 座 头 鲸 游 近 随 机 个 体 去 觅 食 ; 和 否则， 远离 该 
随机 个 体 搜索 其 他 更 适合 猎物 。 

系数 向 量 A 和 C 计算 方法 如 下 : A=2a:r-a、C=2r。 其 
中 ， 随 着 迭代 次 数 的 增加 ，“ 是 从 2 到 0 单调 递减 ， 表 达 式 


pus 


a-2-2L, rAMA, KERAN 


为 [0,1] 。 


bo 包围 收缩 。 座 头 鲸 使 用 种 群 最 优 个 体位 置 来 导航 包围 食 
物 ， 收 缩 捕 食 范围 。 


Xa-X,,-D-A 
D- |c- X -X| 0 


Hp: Xs 是 种 群 中 全 局 最 优 个 体 , 也 是 猎物 所 处 位 置 或 者 
距 猎 物 最 近 的 位 置 ;A、C 同上 ,， 当 | 让 >1 时 ， 表 明 座 头 鲸 远离 随 
机 个 体 Xs ， 向 全 局 最 优 个 体 X us 游 近 。 

cO 螺 线 捕食 。 座 头 鲸 在 向 Xs 游 近 的 同时 , 还 有 一 定 的 概 
率 p 会 以 对 数 螺旋 线 的 轨迹 运动 。 


1+1 (8) 


X, = X, + Dyes xe” xcos(2zl) 
,= |X X] 


其 中 : D, 是 座 头 鲸 当 前 位 置 与 最 优 个 体 的 距离 ;2 为 构造 螺旋 
线 轨迹 的 常数 ;! 是 随机 数 !s[-L1 。 
2.4 混沌 搜索 策略 

在 WOA 算法 中 ， 座 头 鲸 将 最 优 个 体 坐 标 为 导航 ， 并 以 螺 
旋 运 动 同 时 缩小 包围 圈 的 轨迹 游 向 猪 物 。 虽 然 这 可 以 加 快 收敛 
速度 但 是 易 使 群体 快速 聚集 在 一 起 ， 从 而 容易 陷入 局 部 最 优 071。 
混沌 搜索 是 一 种 具有 遍历 性 ， 为 随机 性 特点 的 搜索 方法 。 在 一 
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定 程 度 上 能 够 跳 离 局 部 最 优 值 从 而 找到 全 局 最 优 值 。 但 混沌 搜 
索 无 法 从 经 验 中 学 习 ， 搜 索 不 可 避免 存在 搜索 盲目 性 0 DI, A 
了 结合 WOA 的 快速 寻 优 和 混沌 搜索 的 遍历 性 ， 将 混沌 引入 
WOA 算法 之 中 。 


混沌 搜索 思想 是 : 首先 ， 将 座 头 鲸 的 空间 位 置 通过 混沌 映 
射 转换 为 混沌 变量 ;然后 ,利用 混沌 自 映射 函数 做 遍历 性 搜索 。 


再 将 搜索 得 到 的 混沌 向 量 序列 通过 函数 转换 为 解 向 量 序列 ， 得 
到 最 优 解 向 量 ; 最 后 更 新 座 头 鲸 个 体 空 间 位 置 。 
混沌 搜索 过 程 如 下 : 

a) pL fi s zs [RU or ELS 
D B. H'rpxe(a5). 


四 公式 转换 到 混沌 区 间 〈-1， 


- 2(x-a) | 


Zi 
b; —- a; 


1 (9) 


b) 采 用 比 经 典 罗 辑 映 射 具有 更 好 遍历 性 的 逻辑 自 映射 函数 
20.20 来 产生 多 个 混沌 变量 序列 ， 其 中 Z, ，i=0,1,2,.…; 是 混沌 变 


€, HZeCLD. "42,2005 时 ， 混 沌 就 会 发 生 。 


Z,-1-2Z C10) 


c) 将 所 得 到 的 多 个 混沌 变量 序列 经 下 面 公式 转换 成 多 个 解 
向 量 序列 ， 即 得 到 该 座 头 鲸 个 体 多 个 新 的 空间 位 置 。 


| 1 
x, -—(b —a)xZ,--—(b, +a, 
ios ca) *Z +5% +a) 


d) 对 多 个 空间 位 置 进行 适应 度 值 的 计算 , 选择 最 优 的 适应 
度 值 所 对 应 的 空间 位 置 来 代替 该 座 头 鲸 的 空间 位 置 。 


25 算法 求解 步骤 

a) 参数 初始 化 设置 ， 设 置 种 群 规 模 N， 当 前 迭代 次 数 
最 大 迭代 次 数 ,。。 构 造 螺旋 线 轨 迹 的 常数 b， 随 机 向 量 r， 随 
机 数 1、p。 


b) 初始 化 种 群 ， 按 本 文 提 出 的 初始 化 规则 初始 化 种 群 
计算 每 个 个 体 的 适应 度 值 并 得 到 Xpo o 


c) 如 果 随 机 数 p< zs[0D) 且 |4|<1， 则 个 体 选择 “ 游 走 
砚 食 ”来 更 新 空间 位 置 ， 否 则 个 体 选 择 “包围 收缩 ”来 更 新 空 


间 位 置 。 如 果 随 机 数 p> p* ， 则 个 体 选择 “ 螺 线 捕食 ”来 更 新 空 
间 位 置 ， 并 计算 每 个 个 体 的 适应 度 值 并 更 新 Xs o 


d) 根据 混沌 搜索 策略 ,对 新 更 新 的 群体 中 精英 群体 (适应 
度 值 前 2-10) 进行 混沌 搜索 ， 并 更 新 精英 群体 ， 合 并 未 进行 混 
沌 搜索 的 个 体 ， 从 而 得 到 新 的 群体 并 更 新 Xue 


e) 算法 达到 最 大 迭代 次 数 ， 算 法 终止 ,输出 Xs， 及 对 应 
的 最 佳 调 度 排序 。 


3 ”问题 求解 


3.1 参数 分 析 
WOA 算法 中 


Lu 


FE 要 参数 是 构造 螺旋 线 轨迹 的 常数 pb、 概率 


: vds. 
à soa x 。 ChinaXiv 合 作 期 刊 
REZ, $: 具有 行为 效应 的 含 AGV 有 柔性 车 间 调 度 研 究 


广 、 种 群 数 N、 和 迭代 次 数 。 为 了 研究 b， 广 、N、Tue 取 值 不 同 
对 WOA 的 优化 性 能 的 影响 ， 本 文 对 8x8 测试 示例 [21] 通 过 正 
交 实 验 做 了 以 下 测试 。 
a) 结 合 鲸鱼 算法 相关 文献 参数 取 值 , 确立 参数 水 平 的 取 值 ， 
即 可 建立 四 因素 三 水 平 正 交 实验 表 ( LG) D. 
表 2 参数 水 平 


" 水 平 
参数 1 2 3 
b 0.5 1 2 
p* 0.4 0.5 0.8 
N 100 200 500 
Tmax 50 100 200 


bo 将 各 参数 水 平 组 合 独立 运行 20 次 , 以 20 次 结果 平均 值 
作为 评价 参数 性 能 的 指标 。 其 中 运行 结果 如 表 3 所 示 。 
c) 求 出 每 一 因素 的 极 差 ， 根 据 极 差 ， 得 到 各 参数 的 重要 程 
度 ， 从 而 找到 最 佳 因素 水 平 组 合 。 从 下 表 可 得 出 ，N 的 极 差 最 
大 ， 等 级 最 高 。 说 明 种 群 数目 N 对 算法 寻 优 性 能 影响 最 大 。 然 
而 构造 螺旋 线 轨迹 的 常数 b、 概 率 p 则 对 算法 寻 优 性 能 影响 不 
明显 。 同 样 从 下 图 可 以 很 直观 地 看 出 ， 种 群 数 对 寻 优 性 能 的 影 
响 大 ， 种 群 数 越 大 ， 当 友 代 次 数 相同 时 ， 算 法 获得 的 可 行 解 成 

站 数 增长 ， 因 此 得 到 更 优 解 的 概率 变 大 。 
43 正 交 实验 表 


nin b p* a Tmax es 
1 1 1 1 1 18.55 
2 1 2 2 2 17.8 
3 1 3 3 3 115 
4 2 1 2 3 17.842 
5 2 2 3 1 17.526 
6 2 3 1 2 18.632 
7 3 1 3 2 17.421 
8 3 2 1 3 18.053 
9 3 3 2 1 18.105 
XA 参数 性 能 分 析 
yy 参数 
KF b p* N Tmax 
1 17.833 17.938 18.411 18.061 
2 18 17.793 17.916 17.951 
3 17.86 17.962 17.366 17.682 
极 差 0.167 0.169 1.046 0.379 
等 级 4 3 1 2 


实验 数据 和 分 析 可 得 , 建议 的 参数 组 合 : b=0.5, p' 205, 
N 5500, T, =200 


3.2 算法 寻 优 能 力 验 证 

实验 仿真 环境 为 Windows7，64 操作 系统 ， 处 理 器 主 频 
2.5GHz， 内 存 4G, 采用 MatlabR2014 编译 软件 , 每 个 算 例 运行 
20 次 取 最 优 值 。 群 种 群 规模 NN 为 500 ， 最 大 迭代 次 数 为 200， 
b 取 0.5，p* 取 0.5。 

为 了 验证 粒子 群 算法 的 寻 优 能 力 ， 从 Bilge23 等 提出 的 40 
个 算 例 中 选取 部 分 进行 计算 : 4 台 加 工 设备 、 装 载 台 、 印 载 台 
与 AGV 构成 了 四 种 布局 , 并 且 考 虑 AGV 的 搬运 时 间 和 空 载 时 
间 ， 求 出 最 短 加 工时 间 。 算 例 “Ex21” 表 示 ，2 表示 Jobset2， 
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FPE 


录用 稿 


1 表示 Layoutl. 


加 工 完 成 搬运 至 卸载 台 的 情况 下 ， 采 用 遗传 算法 与 启发 式 混合 
算法 所 得 到 的 结果 。FDE、FMAS 分 别 表示 文献 [25, 26] f th 
法 的 结果 。PSO 和 FPSO 表示 使 用 本 文 改进 粒子 群 算法 ， 分 别 


计算 车 间 调 度 和 和 柔性 车 间 调 度 的 最 短 加 工时 间 。 通 过 实验 结果 


MAS 表示 文献 [23] 提 出 算法 的 结果 。GAHA- 
2. GAHA-3 表示 文献 [24] 提 出 的 2 个 和 3 个 AGV 在 考虑 工件 


China IV 


徐 云 : 具有 行为 效应 的 


8 表明 AGV 的 边际 效用 符合 边际 效 | 


递减 规律 ， 随 着 


AGV 数量 的 增加 ， 总 效用 增加 整个 加 工时 间 减 少 )， 但 边际 


效用 递减 〈 增 加 一 台 AGV 缩短 的 时 间 减 少 )。 


可 知 本 文 提 出 的 改进 粒子 群 人 


法 有 较 好 的 寻 优 能 力 : 运用 本 文 


算法 计算 车 间 调 度 问题 时 ，41.67% 的 解 优 于 MAS (R=2) 与 
GAHA (R22) 的 最 优 解 , 91.67% 的 解 优 于 GAHA (R=3 ) 的 解 。 


在 计算 柔性 车 间 调 度 问 题 


页 时 ,50% 的 解 优 于 FDE(CR=2 ) 与 FIMAS 


(R=2) 的 最 优 解 。 同 时 不 难 发 现 ,增加 AGV 能 够 缩短 整个 完 
工时 间 ， 完 工时 间 越 长 ， 增 加 AGV 更 有 效 。 


€5 算 例 实验 结果 


MAS 
算 例 [22] 
R=2 
EX11 130 
EX12 98 
EX13 109 
EX14 168 
EX21 143 
EX22 86 
EX23 98 
EX24 169 
EX31 142 
EX32 114 
EX33 103 
EX34 167 


FDE FMAS FWOA GAHA WOA 
[24] [25] TX [23] 本 文 


3.3 ”实验 分 析 


3.3.1 AGV 边际 效用 
从 2017-2020 4 


知 ,中 高 端 AGYV 产品 身 
元 。 因 为 AGV 价格 较 高 
因此 本 文 引入 经 济 学 


国 移动 机 器 人 CAGVO. 行业 调研 报告 可 
单价 区 间 在 20-50 万 元 , 甚至 
， 企 业 投资 AGV 数量 决策 尤为 重要 ， 
的 边际 效用 理论 来 衡量 企业 每 投资 一 台 


高 达 100 7j 


AGYV 使 工件 加 


时 间 减 少 的 效果 。 


边际 效用 理 


perps 


为 了 深入 研究 AGV 的 边际 
工 工件 较 多 的 算 例 15 x 10P3)， 


边际 效用 为 4， 而 AGV 数量 
为 边际 效用 小 幅度 增加 的 原 


指出 商品 价值 由 该 商品 的 边际 效用 决定 ， 意 
思 是 每 多 买 一 单位 商品 或 服务 给 消费 者 带 来 的 总 效用 增 量 P。 


根据 理论 ， 本 文 用 多 


em 


所 有 工件 减少 的 时 间 ， 商 品 指 


的 是 AGV， 公 式 如 下 : 


_ ATime 


AU = ATime (12) 


效用 ， 本 文选 取 加 工 设备 和 加 


取 不 同 数量 的 AGV， 分 别 计算 
20 次 ， 所 得 结果 如 下 所 示 。 从 下 表 可 以 发 现 当 AGV 数量 大 于 
等 于 14 时 , 最 优 完工 时 间 在 170。 其 中 , 当 AGV 数量 为 13 时 ， 
EN 14 时 ， 边 际 效用 为 6， 笔者 认 


是 计算 次 数 不 是 足够 多 。 


xET 


X6 基于 不 同 数量 AGV 的 最 优 加 工时 间 

15xi0 1 2 5 6 7 8 5 
PSO[28] 

1 14 
Ald 1555 768 333 282 259 232 2 
WOA 最 优 1325 736 326 287 253 223 205 
时 间 
15x10 10 " 14 15 16 17 18 
PSO[28] 
Sag 199 s 182 184 184 181 18 
WOA 最 优 194 186 170 1 — 170 — 70 — 170 


边际 效用 为 0, 说 明 14 台 和 15 & AGV 总 效用 相同 ， 即 总 效用 


当 AGV-15 时 ， 


在 AGV=14 bi doch 此 时 加 工时 间 主 要 受 加 工 设 备 和 工序 
的 约束 ， 不 能 通过 简单 的 增加 AGYV 来 减少 加 工时 间 。 
800 
589 =a @ a VOA MU 
600 e 
400 z 
00 «6, 
1o 34301811 86460000 
0 99-09. 9.0-0-0-090 0.0 
12345678 9101112131415161718 
图 8 AGV 的 边际 效用 曲线 


3.3.2 学 习 效 应 和 恶化 效应 
D 不 考虑 员工 、 工 件 的 异 质 性 
当 员 工 初始 技能 ， 最 终 技能 相同 且 工 件 开 


始 恶化 时 间 ， 停 


止 恶 化 时 间 相同 时 ， 即 不 考虑 员工 、 工 件 的 异 
同学 习 率 与 不 同 恶 化 因子 对 完工 时 间 的 影响 。 
学 习 曲 线 从 90% 到 60%， 含 AGV 柔性 车 间 调 


质 性 ， 来 讨论 不 
由 下 表 可 看 出 ， 
度 问题 的 完工 时 


间 大 大 缩短 (AGV=2), 即 学 习 率 越 小 , 完工 时 间 越 小 ， 即 学 习 


效应 能 够 提高 工厂 生产 率 。 同 时 可 以 发 现 ， 
加 ， 完 工时 间 逐 渐 增 大 ， 这 表明 较 大 的 恶化 攻 


随 着 恶化 因子 的 增 


子 相当 于 引起 较 


大 学 习 效 果 退 化 。 


表 7 学 习 率 及 恶化 因子 对 寻 优 结果 的 影响 


E31 取 值 最 优 解 
学 习 率 0.9 122 
b-0.1 0.8 111.3004 

0.7 107 
0.6 105.7818 
恶化 因子 0.1 111.3004 
1-0.8 0.15 114 
0.2 119.476 
0.25 122 


20 考虑 员工 、 工 件 的 异 质 性 


i 


以 EX31 [8807 Ui CS 


， 对 所 有 工件 进行 小 批量 生产 ， 


且 批 量 为 1， 2. 3, 4,，5。 下 面 矩 阵 表 示 的 是 每 个 员工 分 别 对 每 


种 工件 的 初始 能 力 、 最 大 能 力 、 学 习 率 ， 如 甜 


行 第 四 列 ) 表示 ， 员 工 2 GTF) 对 工件 2 的 初 


阵 中 0.77 (第 二 
165867] 0.77. KK 


中 工件 的 恶化 因子 为 B=[0.2 02 01 02 015], ， 开 始 恶化 时 
间 t min=[50,50,50,50,50] , W 化 结 X 时 间 


t max-[100,100,100,100,100] r 


0.85 0.78 0.8 0.91 0.78 0.80 0.80 0.70 0.80 0.80 0.70 0.80 
0.92 0.8 0.9 0.77 0.70 0.90 0.90 0.82 0.90 0.90 0.82 0.90 
0.94 0.75 0.7 0.97 0.75 0.70 0.92 0.85 0.70 0.92 0.85 0.70 
0.80 0.6 0.75 0.93 0.80 0.75 0.75 0.70 0.75 0.75 0.70 0.75 


图 9-11 为 当 2 "His Les A OZ, HEITAR A 


位 批量 完工 时 间 的 影响 。 


0.88 0.70 0.80 
0.90 0.82 0.90 
0.92 0.85 0.70 
0.82 0.70 0.75 
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同一 批量 下 ， 只 含 学 


习 效 应 的 单位 批量 完工 时 间 < 


含 学 习 


恶化 效应 的 单位 批量 完 
只 含 恶化 效应 的 单位 批量 完 
工件 开工 时 间 不 长 ， 员 工学 习 还 未 达到 最 


时 间 < 无 效应 的 单位 批量 完工 时 间 < 
时 间 ; 这 说 明 批 量 


z=], AGV=2 时 ， 
大 技能 ， 工 件 处 于 刚 


开始 恶化 阶段 ， 因 此 学 习 效应 起 的 缩短 时 


间 的 效果 之 恶化 效应 


增长 时 间 的 效果 (批量 =2，AGV=2 同 理 )。 

相同 效应 下 ， 批 量 为 2 的 单位 批量 完工 时 间 < 批 量 为 1 的 
单位 批量 完工 时 间 。 这 说 明 , AGV 常 处 于 等 待 工序 加 工 完成 然 
后 去 运输 到 下 一 加 工 设备 的 阶段 , 即 AGV 资源 不 是 瓶颈 资源 。 


其 中 图 


10 可 以 看 出 : 


同一 批量 下 ， 只 含 学 习 


效应 的 单位 批量 完工 时 间 < 无 效应 


的 单位 批量 完工 时 间 < 含 学 习 恶 化 效应 的 单位 批量 完工 时 间 < 


只 含 恶 化 效应 的 单位 批量 完工 
工件 开工 时 间 较 长 ， 员 工学 习 已 达到 最 大 


时 间 ; 这 说 明 批 量 


-3,AGV-2 时 ， 
技能 ， 工 件 处 于 恶化 


阶段 ， 因 此 恶化 效应 增长 时 间 的 效果 > 学 
果 ( 批 量 =3，4，5，AGV=2 同 理 )。 


习 效 应 缩短 时 间 的 效 


不 同 效应 下 ， 批 量 不 同 的 单位 批量 完 


时 间 大 小 关系 在 变 


录用 稿 BEZ, Wi 具有 行为 效应 的 含 AGV 采 性 车 间 调 度 研究 
15. 00 
CN 
13. 00 e^ A E31 r=1 r=2 r=3 
11. 00 * 总 完工 时 间 完工 时 间 / 批 量 总 完工 时 间 完工 时 间 / 批 量 总 完工 时 间 完工 时 间 / 批 量 
批量 =1 186.00 186. 00 112. 00 112. 00 90. 00 90. 00 
ig 09. 00 d 批量 347. 00 173. 50 208. 00 104. 00 151. 00 75. 50 
w 07. 00 PA 批量 =3 507.17 169. 06 304. 00 101. 33 222. 00 14. 00 
E 05. 00 e 批量 =4 680.00 170. 00 390. 00 97.50 295. 00 13.15 
$ 103. 00 e= 批量 -1 批量 -5 864.00 172.80 498. 00 99. 60 365. 00 73. 00 
Ez - =], 
H 101. 00 ACV-2 
IR 99.00 
97. 00 — j| —2 210. 00 
AN Vv No Ve 
ak S sg 2E 170.00 9... ..9..0 
m; 150.00 el 
pe " : P ig 130.00 
图 9 行为 效应 对 单位 批量 完工 时 间 的 影响 H "^ — 9 1=2 
2 110. 
18 $—9 o -9 一 一 -3 
15. 00 SS 
? 70. 00 9——4—4—9 
: ^ d fo 
Ug 110.00 e a aua a 
Š e SESSE 
= 05. 00 erra.. E3, 
H AGV=2 € 11 批量 与 AGV 数量 对 单位 批量 完工 时 间 的 影响 
让 100.00 一 = 。 批 量 =4， 
rae 由 图 11 可 以 看 出 : 
95. 00 — liie 5, 
X^ a em Vis 当 批量 相同 时 ,搬运 机 器 人 越 多 ,单位 批量 完工 时 间 越 短 。 
d b *— LH EH Is L— dm 
x 搬运 机 器 人 为 1 和 2 时 ， 随 着 批量 越 多 ， 单 位 批量 完工 时 间 先 
降低 后 升 高 ， 且 r=1 时 ， 批 量 =3 时 间 最 短 ，r=2 时 ， 批 量 =4 时 
图 10 行为 效应 对 单位 批量 完工 时 间 的 影响 " Eu Lee. . : SEP dro E e 
í me : 间 最 短 。 当 搬运 机 器 人 为 3 时 ， 随 着 批量 越 多 ， 单 位 批量 完 
其 中 由 图 9 可 以 看 出 : 时 间 越 短 。 这 说 明 搬运 机 器 人 为 1 和 2 时 ,由 于 工件 批量 变 多 ， 
AGV 变 成 瓶颈 资源 ， 加 工 设备 需要 等 待 AGV 搬运 工件 才能 


始 加 工 ， 且 搬运 机 器 人 为 1 较 搬 运 机 器 人 为 2 时 ，AGV 变 成 瓶 
搬运 机 器 人 为 3 时 ，AGYV 还 未 变 成 瓶颈 


有 颈 资源 的 批量 更 小 ， 而 


学 习 效应 、 恶 化 效应 现象 如 


客观 存在 ， 同 时 


上 实验 可 知 ， 学 习 效应 、 恶 


生产 制造 过 程 
恶化 效应 在 生产 制 


E 实 际 工 | 


造 过 程 中 对 调度 的 影响 是 不 容 忽 视 的 。 


工序 、 机 床 和 AGV 分 配 的 编码 方式 P1， 用 


4 ”结束 语 
本 文采 用 基于 
3L 维 


向 量 表 示 1 种 可 行 的 生产 调度 。 并 将 混沌 局 部 搜索 策略 运 
用 到 鲸鱼 算法 中 ,有 效 避 免 了 迭代 中 易 陷入 局 


部 最 优 解 的 问题 。 


同时 运用 had 最 后 从 算 例 表 
明了 具有 局 部 混沌 搜索 策略 的 鲸鱼 算法 具有 较 高 的 优化 性 能 ， 


化 。 这 说 明 ，AGYV 资源 约束 和 各 效应 
时 间 。 为 了 剔除 效应 的 影响 ， 进 一 步 探究 
量 之 间 的 内 在 关系 ， 补 充实 验 。 


表 8 批量 与 AGV 数量 对 单位 批 


同 影响 了 单位 批 


完 


AGV 资源 约束 与 批 


完工 时 间 的 影响 


的 因素 还 有 很 多 ， 


发 现 学 习 效 应 能 够 提高 工 | 
学 习 效果 的 退化 ; AGV 资源 约束 和 学 习 效 应 、 
响 了 单位 批量 完工 时 间 ; 随 着 加 工时 1 
工件 加 工时 间 的 效果 越 来 越 明显 ， 直 到 员 
开工 时 间 达 到 工件 恶化 时 间 时 ， 恶 化 效应 延长 工件 
效果 越 来 越 明 显 ， 直 到 超过 最 大 恶化 时 间 。 


比如 遗忘 效应 、 员 工 请 假 、 


生产 M 恶化 效应 反之 ， 相 当 于 


恶化 效应 共同 影 
闻 的 推移 ， 学 习 效应 缩短 
工 能 力 达 到 上 限 ; 当 
加 工时 间 的 
影响 工件 加 工时 间 


ATLAS, AX 


的 角度 分 析 加 工时 


没有 涉及 相关 因素 的 探讨 。 
间 的 影响 因素 。 


姑 此 在 以 后 的 研究 有 必要 从 更 多 


录用 稿 
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